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ABSTRACT 

The properties of galaxies with the reddest observed R — K colors (Extremely Red Objects, 
EROs), including their apparent division into passive and obscured active objects with roughly 
similar number densities, are a known challenge for models of galaxy formation. In this paper 
we produce mock catalogues generated by interfacing the predictions of the semi-analytical 
MORGANA model for the evolution of galaxies in a ACDM cosmology with the spectro- 
photometric + radiative transfer code GRASIL and Infrared (IR) template library to show that 
the model correctly reproduces number counts, redshift distributions and active fractions of 
R — K > 5 sources. We test the robustness of our results against different dust attenuations 
and, most importantly, against the inclusion of TP-AGB stars in Simple Stellar Populations 
(SSPs) used to generate galaxy spectra, and find that the inclusion of TP-AGBs has a relevant 
effect, in that it allows to increase by a large factor the number of very red active objects at 
all color cuts. We find that though the most passive and the most obscured active galaxies 
have a higher probability of being selected as EROs, many EROs have intermediate proper- 
ties and the population does not show bimodality in specific star formation rate (SSFR). We 
predict that deep observations in the Far-IR, from 100 to 500 fim, are the most efficient way 
to constrain the SSFR of these objects; we give predictions for future Herschel observations, 
and show that a few objects will be detected in deep fields at best. Finally, we test whether a 
simple evolutionary sequence for the formation of z = massive galaxies, going through a 
sub-mm-bright phase and then a ERO phase, are typical in this galaxy formation model. We 
find that this sequence holds for ~ 25 per cent of z = massive galaxies, while the model 
typically shows a more complex connection between sub-mm, ERO and massive galaxies. 
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1 INTRODUCTION 

The wealth of data coming from recent multiwavelength surveys 
has made it possible to study in detail distant galaxy populations 
from their UV-optical-Infrared broad-band spectral energy distri- 
butions (SED). Among the various categories that have been obser- 
vationally defined by means of color selection techniques, the Ex- 
tremely Red Objects (EROs) have attracted wide interest. Defined 
on the basis of of their very red R — K color (in the Vega system 
R—K > 5 at least), th ey were first addres s ed in the context of dee p 
near-infrared surveys felston et alj|l988l : iMcCarthv et all 1 19921) . 
Such red colors are expected to arise as a combination of an intrinsi- 
cally red SED and a large K-correction at 1 < z < 2. Then, EROS 
wer e regarded as natural counterparts of local elliptical galaxies 
(e.g. lCimatti et alj|2002af) . Moreover, their relatively large /\-band 
fluxes suggested high stellar masses, comparable with local ellip- 
ticals, leading to the interpretation of EROs as already assembled 
elliptical galaxies and the redshift range z > 2 as their formation 
epoch. The correct prediction of such an early and efficient forma- 



tion of massive galaxies at high redshift is a well known challenge 
for models of galaxy formation based on the CDM cosmology: in 
fact early implemen tations predicted lower space densit ies of EROs 
than observed fe.g jFirth et alj2002l ; ISmith et alj|2002l) . 

On the other hand, it soon became evident that the same red 
R — K colors might be linked to a completely different class 
of objects such as the active (starburst) galaxies embedded in 
dust shells that cause high extinction of the SEDs JCimatti et al.l 
l2002al : |Pozzetti & Mannuccil2000l : ls"mail et alj2002l) . Interestingly 
enough, the relative contribution to the total space density by the 
two sub-popul ations is similar, w ith a passive fraction of 48 per cent 

at K < 19.2 l lKong et al. 2009, but a ratio between 30 and 70 per 
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cent has been reported also by Cimatti et al. 2002a and Smail et al. 

l2002h , and 58 per cent at K < 20.3 jMivazaki et"ai]|2003h . How- 
ever, splitting the observed ERO population in two is a difficult 
task. The cleanest way relies on high signal-to-noise spectroscopy, 
able to reveal either emission lines linked to ongoing star forma- 
tion or the spectral breaks typical of aged stellar populations (see 
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e.g JCimatti et al .l2002ah . Alternatively, the presence of a starburst 
may be revealed by additional im aging in the far-in frared, where 
the absorbed light is reprocessed dKong et al.ll2009h, or indirectly 
through the morph ology of these objects I Stiavelli & TreulfeoOll ; 
ICimatti et al]|2003l) . Presently, all these techniques can be success- 
fully applied only to a small fraction of the overall population, typ- 
ically limited to the brightest sources. 

As mentioned above, the redshift distribution and the expected 
physical properties of the two ERO sub-populations have been a 
long-standing problem for theoretical mode ls of galaxy formation 
dCimatti et al.ll2002al ; [Somerville et alj2004) . More recent models, 
based on improved treatment of baryonic physics (|De Lucia et al.l 
2004 iBowei- et al.l 120061 : iMenci et all |2006t iMonaco et al.ll2007F 



Lagos et al.ll2008l) " are able to give a better representation of the 



formation and assembly of massive galaxies at z < 3. Only a few 
of these models have been directly tested by comparing the pre- 
dicted properties of EROs with data: in most cases a reasonable 
success in reproducing the abundance of passive EROs has been 
claimed, giving support to the interpretation of these objects as 
passive massive galaxies at z ~ 1 — 2 that are the progenitors 
of local elliptical galaxy. At the same time, the same models do 
not produce a sufficien t number of active EROs. More specifically, 
iNagamine et al.l (2005) implemented sub-grid baryonic physics in 
both a Lagrangian and Eulerian code; they found the correct abun- 
dance of massive galaxies at z ~ 2, but could find almost no 
passive ERO, and in order to produce enough red galaxies they 
were forced to assume a dust attenuation as high as an equiva- 
lent E (B — V) = 0.4 for the whole mock catalogue. iKang et al .1 
J2006h correctly reproduced the space density of massive galaxies 
at z ~ 2, but underpredicte d the number of ver y red galaxies, in 
particular of the active class. IMenci et al.l J2006) proposed a more 
successful model able to reproduce the color distribution of galax- 
ies at z < 1 and the space density of massive galaxies at 1 < z < 2 
at the same time. This model is able to reproduce the ERO num- 
ber counts and redshift distribution, but it predicts a relatively high 
fractio n of passive objects (~80 p er cent at K < 22). More re- 
cently, [Gonzaifiz^Perezet^^ implemen- 
tations of the GALFORM model l lCole et alj|2000T) to observational 
data o f the ERO population. Their best fit model (i.e jBower et al .1 
120061) produced the correct number of EROs, even when a very red 
cut (up to R — K > 7) was used, but again the predicted ERO popu- 
lation was dominated by passive objects, with only a small fractio n 
of active galaxies. On the other hand the model by iBaughl J2006I) . 
which gives a good fit of the sub-mm galaxy population, was found 
to significantly underestimate ERO number counts. 

A remarkable attempt to put quasars, bright sub-mm galaxies, 
E ROs, and low-z ellipti cals into a unified scenario was proposed 
by iGranato et al. I l l2004l) . In their model a typical low-z elliptical 
forms the bulk of its mass in an intense obscured starburst phase at 
2 ~ 2, when it is visible in the sub-mm band. In the galaxy nucleus, 
star formation triggers accretion onto a seed black hole; the end of 
the star formation phase is caused by feedback from the resulting 
quasar, able to wipe away the interstellar medium from the galaxy 
and quench star formation. After ~ 1 Gyr the galaxy gets very red 
and becomes an ERO. It then evolves passively into a typical local 
elliptical galaxy. 

The above mentioned papers have clarified the role of the ERO 
population as a powerful test for theoretical models of galaxy for- 
mation. These should reproduce both their global number and the 
diversity of their physical properties. A successful model would be 
helpful in investigating the relation of EROs with other observa- 
tionally selected galaxy populations, in particular with the sub-mm 



galaxies. Apparently, the most challenging requirement for pub- 
lished model is to produce a high fraction of active EROs, espe- 
cially when deep samples are concerned. In principle, active EROs 
should be easy to produce in a hierarchical model of galaxy for- 
mation, because the redshift range where we observe them cor- 
responds to the p eak of the cosmic star for mation rate density. 
iKong et al.l d2006l) and lNagamine et alj l l2005h interpret the tension 
between their models and data as a sign of over-simplified recipes 
in treatin g baryonic physics or d ust attenuations. Regarding the lat 
ter point. 
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for dust attenuation, calibrated using low-z samples, provide very 
scattered extinction values at z ~ 2 and, more interestingly, do 
not match the predictions of the GRASIL code dSilva et aTll 19981) 
that explicitly solves the equations of radiative transfer. Another 
important point is the contribution of th e TP-AGB stars to the syn- 
thetic SEDs of model galaxies (see e.g. iMar aston 2005), which is 
not included in the a bove mentioned models. It has been shown by 
iTonini et al .1 d2009al lbT) that the inclusion of this extreme phase of 
stellar evolution has a strong impact on the expected photometry 
and colors of simulated galaxies at 1 < z < 4, where intermediate- 
age stellar populations dominate the restframe A'-band starlight. In 
the case of EROs, an increase of their number density is expected 
when TP-AGB s are included due to the reddening of rest-frame 
V — K. All these elements critically affect the modeling of galactic 
SEDs: they are therefore a primary concern for the interpretation of 

the physical ingredients of models. 

In this paper we use the MORGANA model dMonaco et al.l 

as updated in iLo Faro et all l2009h , interfaced with the 
GRAS IL spectro-photometric + radiative trans fer code Silva et al.l 
dl998l) and to Infrared (IR) template library dRieke et alj|2009i) . 
to produce mock s amples of EROs. We bo th use the Padova 
dBertellietal1ll994h and the iMarastor] d2005l. M05 hereafter) li- 
brary of SSPs. We find reasonable agreement with data in terms of 
number counts, redshift distribution and fraction of active objects, 
especially when M05 library is used (details on the effect of this 
modification on the whole galaxy population will be the subject of 
a forthcoming paper). This prompts us to use the model to inves- 
tigate the physical properties of these object. We address whether 
there is a true active/passive bimodality of the ERO population and 
what is the connection between EROs, sub-mm galaxies and mas- 
sive galaxies at z — 0. We finally give predictions of ERO proper- 
ties for future surveys with the Herschel satellite. 

This paper is organized as follows. In sec.|2]we recall the main 
feature of the theoretical models MORGANA and GRASIL. In sec. [3] 
we present our results: in sec. |3.1| we analyze the properties of the 
model ERO sample and in sec. 13.21 we discuss the connection be- 
tween the ERO population, the sub-mm bright galaxies and the 
massive galaxies at z = 0. Finally in sec. [4] we give our conclu- 
sions. Throughout the paper we assume magnitudes are in the Vega 
system (unless otherwise stated), and a cosmological model con- 
sistent with WMAP3 r esults (Q = 24, Q A = 0.76, h = 0.72, 
as = 0.8, n sp = 0.96. lKomatsu et alj|2009h . 



2 MODEL 

We use the semi-analytical MOdel for the Rise of GAlax i es aN d 
AGNs (MORG ANA), originally pre sented in lMonaco et al. I d2007h . 
and updated in lLo Faro et alj d2009l) . where it was optimized for the 
WMA P3 ACDM cosm ology used in this paper, Padova SSP library 
and a lChabrieJ d2003l) stellar Initial Mass Function (IMF). All de- 
tails are reported in the former paper, while a brief outline is given 
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in the latter; we refer to them for a description of the model. The 
only new extension we present in this paper is the use of M05 li- 
brary as an alternative to the Padova ones. The predictions of MOR- 
GANA we re compared with a number of obs ervations in a series 
of papers dFontanot etal]|2006ll2007ll2009bh . In those papers the 
model was calibrated to reproduce the local stellar mass function 
and the star formation rate density as a function of redshift, while 
some parameters (especially those regulating accretion onto central 
black holes) were fixed by requiring a good fit of AGN luminos- 
ity functions. Among the most interesting successes of the model, 
MORGANA is able to correctly reproduce, with the same combina- 
tion of parameters, (i) the evolution of the AGN population and in 
particular the optical and X-ray luminosity functions; (ii) the red- 
shift distribution and luminosity functions of A-selected samples 
in the redshift range < z < 3; (iii) the number counts of 850^im- 
selected sources assuming a conservative Salpeter or Chabrier IMF. 
Moreover, the progressive build-up of bright Lyman-break galaxies 
at z > 4 is recovered after the optimization of dust attenuation pa- 
rameters that have a dramatic influence on rest-frame UV radiation. 

jf-ba nd and sub-mm coun ts of MORGANA galaxies were 
shown by iFontanot et al.l |2007) for the original version of the 
model. With the change of c osmology, IMF and the improvements 
given in iLo Faro et aflfc009l) . A -band number co unts are very sim- 



ilar to those presented in Fontanot et al.l d2009bh . with a better fit 
at K ~ 18, roughly corresponding to the knee of the local lu- 
minosity function. The agreement of the redshift distribution of 
K < 20 galaxies with observations shows a slight but appreciable 
improvement. Though the stellar mass fun ction is well recovere d 
at high masses with the improved model dFontanot et alj|2009bl) . 
we still produce an excess of bright galaxies in the A-band local 
luminosity function (for K < —25.5, amounting to ~ 1 dex at 
Mk = —26.5). This is due to the too high specific star formation 
rate of local massive galaxies (see below), which results in mass- 
to-light ratios M±/Lk ~ 0.5 Mp) /La, to be co mpared, e.g., with 
the value of 0.75 assumed bv lCole et al.l d200ll) to convert the lo- 
cal A-band luminos ity function into a stellar mass function (see 
IFontanot et alj|2007l for a complete discussion of this issue). Sub- 
mm counts are depressed by a factor of 0.5 dex when passing from 
Salpeter to Chabrier IMF, while their redshift distribution does not 
change appreciably: this implies that in the new model sub-mm 
number counts are consistent with data up to a flux of 3 mJy, while 
they are undepredicted at brighter fluxes. 

Previous papers have found a number of points of severe ten- 
sion between model and data, some of which are also present for 
most similar models. 

(i) The colors of mod el massive galaxies at low redshift are too 
blue dKimm et alj|2009l) : this is a peculiar problem of MORGANA 
and is mainly due to the modeling of AGN feedback: because star 
formation is required to trigger black hole accretion, an efficient 
AGN feedback requires some small amount of star formation to 
take place, and this makes massive galaxies bluer than observed. 

(ii) The color distribution of low-z satellite galaxies is not cor- 
rectly reproduced, with the frac tion of red satellites ou tnumbering 
the corresponding blue fraction dWeinmann et al.l200db . This is di- 
rectly related to "strangulation", i.e. the assumption that satellites 
are stripped of all t heir hot gas when they infall on a larger group. 
iKimm et all d2009l) recently showed that this problem is common 
to most semi-analytic models of galaxy formation. 

(iii) The so-called "downsizing" trend in galaxy formation (e.g. 
more massive galaxies forming on a shorter timescale and at higher 
redshift with respect to lower mass counterparts) is not fully re- 



covered. The tension between model predictions and observations 
is significant for low-to-intermediate stellar mass ga laxies (M« < 
10 11 M a): they form too efficiently at high redshift dLo Faro et al.l 
120091 see also) and are too passive at later time, so that the model 
overpredicts their observed stellar mass function at z > 1. As a 
result such small galaxies are predicted to host too old stellar pop- 
ulation at z = 0. IFontanot et al. I d2009bl) showed that this behav- 
ior is common to most published mode ls based on ACDM (i.e. 
IWang et alj|2008l ; ISomerville et alj200l) . 

(iv) The number of massive galaxies at z > 3 is underpre- 
dicted by MORGANA: despite the evolution of the A'-band lumi- 
nosity function is rea s onably reproduced by the mo del up to z ~ 2 
dFontanot et al ]|2007h . lCirasuolo et al.ld2007ll2008h showed that at 
higher redshifts the bright end is systematically underpredicted. A 
similar effect is seen i n the evolution of the stellar mass function 
dFontanot et al. I l2009bl) . This tension is expected to weaken when 
TP-AGB stars are inc luded in the popula tion synt hesis model used 
to gen erate the SEDs dTonini et alj2009bh . fhough lMarchesini et al.l 
d2009l) find that this is not sufficient to remove the discrepancy. 

(v) Specific star formation rates of active galaxies at z ~ 2 in 
MORGANA galaxies are higher than those of other mo dels, but they 
may still be a factor of ~ 3 l ower than observed ones dSantini et al .1 
120091 : IFontanot et al.ll2009bh . though this difference may well be 
due to systematics in star formation rate estimators. 

(vi) Higher specific star formation rates are clearly the key 
for the success of MORGANA in reproducing the elusive sub-mm 
counts, but the redshift distribution of the most lu minous sub-mm 
source s peaks at lower redshift than observed by IChapman et al] 
d2005h . This discrepancy is clearly connected with the lack of 
downsizing in the model: massive starbursts should take place ear- 
lier than average, not later as in our model. 

All these issues can influence our results on the ERO popula- 
tion; for instance, a population of strangulated satellites could arti- 
ficially boost the number of EROs, while residual star formation in 
massive galaxies could depress it. We will keep in mind all these 
discrepancies when interpreting our results. 



2.1 Generation of mock catalogues 

We generate our m ain mock catalogue using a PINOCCHIO 
dMonaco et al.l200"2h simulation of a 200Afpc comoving box, sam- 
pled with 1000 3 particles. The particle mass is then ~ 1O 9 M0 and 
the smallest halo we consider contains 50 particles, for a mass of 
5.1 x 10 10 Mq. The typical stellar mass of the central galaxy con- 
tained in the smallest DM halo at 2 = is ~ 5 x 1Q 8 M(7). The 
same s ampling strategy to generate deep fields as in lFontanot et al .1 
d2007l) is adopted here, with the only difference that, in order to 
have a sufficient number of EROs, all merger trees in the box are 
used (wtvcc = 1 in the original notation, see the original paper for 
more details on sampling strategy and catalogue generation). The 
final reference mock catalogue consists in ~ 240000 model galax- 
ies, covering a redshift range < z < 8. 

MORGANA provides predictions on galaxies in terms of phys- 
ical, non observable quantities. We then interface MORGANA 
with the spectro- photometric + radiative transfer code GRASIL 
dSilva et al.|[l998h . This code computes the synthetic UV-to-radio 
SEDs of galaxies from the known star formation history, cold gas 
mass, cold gas metal mass and size, separately provided for a bulge 
and a disc component. Further assumptions are made on bulge pro- 
file, vertical width of the disc, dust distribution in molecular clouds 
and diffuse gas. Stars are assumed to be born within the optically 
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Figu re 1. ERO cumulati ve number counts in the reference mode l for four R— K (Veg a) color cuts , as giv en in the panels. Observat ional datapoints are taken 
from lDaddi et alj |2OO0l blue empty circ les), Ismail et alj j2002[ cyan empty triangles). ISm ith et al. J2002L ma genta filled triangles ), iRoche et all j20O3l black 
empty diamonds). iMiyazaki et alj j2003l magenta asterisks), yiusanenjfe_J^anssoril d2004l red filled squ ares),| Brown e t alj j2005l black stars), ISimpson et alj 
, KormeTaTl 120061 green filled circles), LawrencTet*alTl2*007r cvan empty squares). [Conselice et al.l|2008t black filled diamonds). Solid 



12006, green crosses), 

black lines refer to the prediction for the whole ERO population, red dashed and blue dot-dashed lines refer to the passive and active classes respectively. 



thick molecular clouds and to gradually escape from them as they 
get older. This gives rise to an age-selective extinction, with the 
youngest and most luminous stars suffering larger dust extinction 
than older ones. The radiative transfer of starlight through dust is 
computed assuming an axisymmetric bulge-disc system, yielding 
the emerging SED. Some of the parameters needed by GRASIL 
are not directly provi ded by MORGANA. We fix their values as in 
iFontanot et alj 120071) . where a choice was made which is good for 



the average galaxy population: in particular we use an escape time- 
scale of young stars for the parent molecular clouds equal to 10 7 yr 
and a 0.5 fraction of cold gas mass in the for m of molecular clouds . 
These values differ from the choice done in iLo Faro et all ( l2009t) . 
which was optimized for high-z Lyman-break galaxies: escape time 
equal to 0.3 x 10 7 yr and a 0.95 molecular fraction. As shown by 
ISilva et alj J200I) . the final fraction of active EROs in their model 
depends strongly on these unconstrained parameters: this is espe- 
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Figure 2. ERO cumulative number counts in the model using M05 SSP library for four R— K (Vega) color cuts. Symbols and lines are as in figure[T] 



cially true for the objects at z > 1, which require a detailed com- 
putation of the highly uncertain dust obscuration in the restframe 
UV. They reported that obscuration is especially sensitive to the as- 
sumed molecular fraction. Because our aim in this paper is to test 
the predictions of already published versions of the model and not 
to find the best-fit parameters for this specific obser vables, we de- 
cided to use the parameters of iFontanot et alj d2007l), and test how 
results change using the high molecular fraction of ILo Faro et al.1 
d2009h . Apparent magnitudes are computed by adding a fiducial 
0.1 mag random Gaussian error. Nu mber counts and redshi ft dis- 
tributions are then constructed as in IFontanot et alj {2007). It is 
worth stressing that our reference MORGANA+GRASIL model is 



calibrated using the Padova SSP library. Here, we a lso produce a 
catalo gue using the M05 SSP library (which adopts a lKroupa et al.1 
i 19931) IMF - we neglect the small differences with the Chabrier 
IMF), to check the effect of the inclusion of the TP-AGBs on the 
ERO statistics. 

Given the low space density of the ERO population, a higher 
sampling of model gal axies is required with respect to the mock 



catalogues defined in IFontanot et al 



up computations (see IFontanot et al 



d2007t) . In order to speed 
120071 for a general discus- 



sion about the time requirements of the MORGANA+GRASIL al- 
gorithm), we modify our standard approach as follows. For each 
model galaxy we compute the UV-to-Near-IR SED with GRASIL, 
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including dust attenuation but excluding the computation of the 
Far-IR emission due to heated dust, which is by far the most time- 
consuming step. With this choice we are then able to generate the 
synthetic SEDs in a few seconds instead of ~ 5 — 10 min. We 
then compute the total absorbed starli ght and we use it to select 
an appropriate IR template from the iRieke et all l l2009l) library. 
We renormalize the chosen template to the total absorbed starlight 
and use it to model IR emission, which is sig nificant for fluxes at 
A > 3/im. As shown in lFontanot et al.1 d201(j|) . this procedure pro- 
vides an estimate for the SED which is however, for high-redshift 
galaxies, biased low by ~ 0.2 dex at high sub-mm luminosities but 
is correct on average at fluxes of < 1 mjy, with respect to GRASIL 
code. As a cross check, we apply a full GRASIL calculation to a 
much smaller mock catalogue of ~ 32000 galaxies and compare 
the resulting number counts and redshift distribution with those 
based on IR templates. No significant differences were found at 
A < 850^im. The normalization of sub-mm number counts is cor- 
rectly reproduced at fluxes /sso/jm ~ 0.5 mjy: at brighter fluxes 
the counts are underpredicted by a factor of 4 at /sso/jm ~ 3 mjy, 
due to the low biasing of this procedure with respect to the predic- 
tions of the full GRASIL calculations, while at fainter fluxes (> 0.1 
mjy) they are overpredicted by a similar amount. 



3 RESULTS 



In fig.Q]we compare the predicted ERO cumulative number counts 



(black line) with available data (Daddi et al. 2000; Smith et al. 


2002 


;ISmail et alJl2002l;lRoche et alj|2003; 


Vaisanen & Johansson 


2004 


; iBrown et alj|2005l; iKons et all 20061: ISinroson et alj|2006l: 


Lawrence et al. 2007: Conselice et alj 2008b. for color cuts (R — 



K) > 5.0, 5.3, 6.0 and 7.0. We first consider the predictions of 



our reference model: for the first color cut, the agreement is good 
down to K ~ 20, where the model starts to overpredict number 
counts. This behavior is very similar to that of the total sample (fig- 
ure 2 of lFontanot et alj2 007). We stress that this level of agreement 
is obtained, with no para met er tuning, from the s ame model used 
in lFontanot et al.ld2009bT) and lLo Faro et al.l ( l200Sl) . The agreement 
worsen at redder colors: the brightest sources at (R — K) > 5.3 are 
underpredicted by an order of magnitude, lowering to a factor of 3 
at K ~ 18; we find an overall underestimate by an order of magni- 
tude for R — K > 6 and a much more severe one for R — K > 7. 

We then split the whole ERO population into active and pas- 
sive subsamples using the instantaneous specific star formation, 
SSFR = SFR/M*. As a threshold between the two popula- 
tions we assume the same criterion as in lBrinchmann et akl \ 20041 
SSFR = 10" 11 2/r _1 ) at all redshifts: at the depth of K = 22 and 
for R — K > 5.0 we classify ~ 68 and ~ 32 per cent of model 
EROs as active and passive respectively. At K < 20.3 the relative 
fractions of active and passive EROs b ecome ~ 62 and ~ 38 per 
cent to be compared with the results of lMivazaki etalj d2003l . , 42 
and 58 per cent respectively), while at K < 19.2 the fract i ons are 
~ 53 and ~ 47, consistent with the results o f lCimattietaII ( l2002al 
both 50 ± 20 per cent). Finally at R - K > 5.3 and K < 20.5 our 
model predicts a relative fraction of 56% ac tive EROs, to be com- 
pared with the results of ISmail et alj d2002l) . who found an active 
fraction ranging from a lower limits of 30% up to an upper limit 
of 60%. Figure [T]gives the predicted cumulative number counts for 
the two classes of EROs; for the first color cut passive EROs (red 
dashed line in the figure) are more numerous at bright fluxes, active 
ones (blue dot-dashed line) start to dominate at K > 18.5. 

To test the stability of our results against different SSP li- 
braries and models of dust attenuation we repeat the analysis using 



different assumptions for these ingredients. We create three mock 
catalogues by interfacing the same MORGANA realization with the 
following modifications to our standard GRASIL setup: 

(i) we consider a higher molecular fraction (0.95) and a lower 
escape time for sta rs (3 Myr) in active galaxies, as used by 
iLo Faro et ail d2009h to reproduce the properties of bright Lyman 
Break Galaxies; 

(ii) we model dust extin ction using the prescription proposed by 
iDe Lucia & Blaizo tl ( 120071) . which combines their fitting formula 
for the optical depth in the V-band with a slab geometry and a com- 
posite age-dependent attenuation law (i.e. assumin g a Milky Way 
extinc tion curve for the diffuse "cirrus" dust and the lCharlot & Falll 
|2000| power-law attenuation law for the younger stars). 

(iii) we use the M05 SSP library in place of the Padova one. 

Our results show that changing dust parameters or the modeling of 
extinction has only a marginal effect on the shape and normaliza- 
tion of the number counts. Increasing the fraction of cold gas in the 
molecular phase depresses the counts of active objects by 0.3 dex 
at K ~ 20, while the passive population is unaffected. A similar 
decrease ( 0.1 dex) is observed when using the analytical attenua- 
tion law of IDe Lucia & Blaizotl d2007h : in this case only ~ 49 per 
cent of the original EROs are selected, while an additional ~ 21 per 
cent (with respect to the original number) is added to the sample. 
In both tests on dust attenuation, all the removed objects, as well 
as the new inclusions, belong to the active subsample. Since the in- 
trinsic synthetic SEDs do not change in these tests, the differences 
in the colors of model galaxies aris e from the different att enuation 
laws. Indeed it has been shown bv lFontanot et alj d2009al) that De 
Lucia & Blaizot's analytical attenuation law is a good approxima- 
tion of the mean attenuation law in a MORGANA+GRASIL realiza- 
tion, but the scatter is large. At the same time, a different choice of 
dust parameters is able to change both the mean and scatter of the 
distribution. 

On the other hand, we find that using M05 SSP library gives 
the strongest effect on ERO number counts, causing a relevant in- 
crease of EROs at all color cuts; the active population becomes 
dominant at all fluxes. In this case the change in the colors of model 
galaxies is mainly driven by a change in the intrinsic shape of their 
synthetic SEDs. In fig.[2]we show the integral number counts in this 
last model: R — K > 5.0 counts are slightly overpredicted at all 
magnitudes, while the 7? — K > 5.3 bright counts now agree with 
observations; more interestingly the R — K > 6.0 are now in excel- 
lent agreement with observations, while the R — K > 7.0 sources 
are still under represented in the model. Using a higher molecu- 
lar fraction or De Lucia & Blaizot's analytical attenuation law in 
conjunction with the M05 SSP library allows a further (marginal) 
improvement of model predictions with data. 

In fig. [3] we compare the redshift distribution of model 
EROs with available data. In the left panels we consider the red- 
shift distribution of t he global populati on dCimatti et alj l2002hl ; 
iMivazaki et all 120031 : ICaputi et alj|2004l) different Tf-band mag- 
nitude limits. For our standard model (solid line), the agreement 
is good in the probed range: EROs are mostly present in the red- 
shift range [0.8, 2] (mean redshift 1.34 at K < 20), with a sharp 
low-redshift cutoff which is well reproduced and a high-redshift 
tail which is under-represented by the model. Regarding the nor- 
malization, we notice again some excess at fainter magnitudes and 
some possible underestimate at K < 19.2, which is however not 
confirmed by other surveys shown in figure Q] When the sample is 
split into active and passive populations, we see a marginal trend 
of passive galaxies peaking at slightly lower redshift, which our 
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Figure 4. Upper panel: 2D probability density of model galaxies in the 
Mi, - SSFR space. Blue continuous lines are iso-density contours spaced 
by factors of 10 in density, and are relative to the control sample of the all 
galaxies with K < 22 and in the redshift range z = [1,2]. Dashed red lines 
give the same quantity for the ERO population. Gray-scale contours give the 
fraction of EROs with respect to the control sample in the same Af*-SSFR 
bin, using the scale given on the right. The mid and lower panels report the 
same quantity for central and satellite galaxies respectively. 
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Figure 5. 2D probability density of model galaxies in the SSFR- E(B—V) 
space. Symbols are as in fig. [4] 



All Galaxies 



0.5 r 



0.4 r 




I 



0.90 _ 
0.85 _ 
0.80 _ 
0.7S _ 
0.70 _ 
0.65 _ 
0.60 _ 
0.55 _ 
0.50 _ 
0.45 _ 
0.40 _ 
0.35 _ 
0.30 _ 
0.25 _ 
0.20 _ 
0.15 _ 
0.10 _ 



0.0 0.2 0.4 0.6 0.8 1.0 
Age galaxy/Age Universe 



1.2 



Figure 6. 2D probability density of model galaxies in the space defined by 
mass weighted age, divided by the age of the Universe at the redshift of 
observation, and E(B — V). Symbols are as in fig. [4] 



model roughly reproduces (mean redshift 1.27 and 1.13 at a depth 
of K = 20 for the acti ve and passive popula tions respectively); 
moreover, from the riche rlMivazaki et aljj2003h sample we see that 
the under-reproduction of the high-redshift tail is mostly due to pas- 
sive galaxies, while active EROs are possibly over-represented. The 
modified models with different dust parameters or analytical atten- 
uation law predicts redshift distributions that are almost undistin- 
guishable from the standard model. The model with the M05 SSPs 
predicts similar redshift distribution (mean redshift 1.34, 1.25, 1.15 
at K < 20 for the total, active and passive populations respec- 
tively), but higher number densities (dashed lines). This is espe- 
cially relevant for the active subpopulation. We als o show the red- 
shift distribution for the redder cut R — K > 5.3 dSimpson et al.l 
120061) and we find that they preferentially select EROs at higher 
redshift; our model with M05 library reproduces this trend. 

These tests show clearly that, among the three changes in SED 
modeling that we tested, the inclusion of TP-AGBs has a strongest 
effect on the number counts, and in particular, on the active sub- 
population. The detailed analysis of the effects of this inclusion on 
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Figure 7. Upper panel: 2D probability density of model galaxies in the M* 
- Z+ space. Symbols are as in fig. [4] The mid and lower panels report the 
same quantity for passive and active galaxies respectively. 



the whole galaxy population and its redshift evolution is clearly be- 
yond the aim of this paper and we leave it to future work. For this 
reason, and in order to allow a direct comparison with the previ- 
ously published results for the same model, we decided to discuss 
the physical properties of EROs only for the reference model with 
Padova SSPs, which is able to reasonably reproduce counts and 
redshift distributions of the R — K > 5.0 sources, even when we 
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split the population into active and passive subsamples. We checked 
that the results presented in the rest of this paper do not change sig- 
nificantly when the model including M05 SSPs is used. 

Summing up, the three main points of disagreement of our 
reference model with data are: (i) an excess of faint EROs, (ii) an 
underestimate of EROs at high color cuts, (iii) an underestimate of 
high-redshift passive EROs. These discrepancies are consequences 
of the known problems of the model: the excess of small galax- 
ies at z ~ 1 that affects most galaxy formation models reflects 
in the excess of faint EROs, while because the reddest galaxies 
are higher-redshift sources, points (ii) and (iii) are usually related 
with the dearth of h igh-redshift massive galaxies in MORGANA 
jFontanot et al .120071) . It is worth stressing that when TP-AGB stars 
are taken into account, the tension at points (ii) and (iii) is consid- 
erably reduced (but not completely solved): this is an hint of the 
relevance of the detailed inclusion of this stellar evo lutionary phase 
in the algorithms for the estimate of stellar masses iMaraston et al .1 
l200d) . 

3.1 Properties of the ERO populations. 

The good level of agreement of model EROs with observations at 
R — K > 5, especially in terms of fraction of active and pas- 
sive objects, encourages us to use the model to get insight on the 
properties of the two ERO populations and on their relation with 
sub-mm galaxies and local massive galaxies. To this aim we con- 
sider the predicted physical properties of the K < 22, R — K > 5, 
1 < z < 2 ERO galaxies, and compare them with a control sam- 
ple containing all galaxies with K < 22 and in the same redshift 
range. In the following we will consider only predictions relative to 
our reference model, but we checked that the inclusion of any of the 
modifications we discussed in the previous section does not affect 
our conclusions. In figure|4]we show the 2D probability density of 
ERO (red dashed lines) and control sources (blue continuous lines) 
in the stellar mass - SSFR space. The gray-scale contours give the 
fraction of galaxies in that stellar mass and SSFR bin that meet the 
selection criterion. Overall, EROs do not populate a special region 
in this space with respect t o the global populatio n at sim ilar red- 
shift and, consistently with I Gonzalez-Perez et al.l J2009I) . they do 
not show a bimodal distribution of either passive or active galaxies. 
Passive galaxies show a higher probability of satisfying the ERO 
selection; the most active galaxies are typically selected out. 

A crucial point of galaxy formation models is the different be- 
havior of central and satellite galaxies; in particular, strangulation 
of satellites may artificially boost the number of EROs. To deepen 
this point we separately show the ERO and control distributions of 
centrals and satellites in the M + - SSFR space, in the two upper 
panels of fig. [4] As expected, satellites are more passive than cen- 
trals and centrals are more massive than satellites, and this is true 
both for the ERO and for the control samples, but the two popula- 
tions of active and passive EROs are present in both sub-samples, 
and the differences between the distributions in the two cases are 
not so strong. We cross-checked this result by comparing integral 
number counts of ERO central galaxies with all EROs, finding only 
a small difference. In particular, the faint end of the number counts 
is reduced by < 0.3 dex, while the bright end is unaffected. This 
is due to the fact that the population of M* > 10 10 Mq galaxies 
is dominated by centrals. We conclude that the statistics of model 
EROs is unaffected by artificial strangulation of satellites. 

The absence of a clear bimodality in the distribution of ERO 
SSFR is at variance with early ideas on the nature of these objects 
as either massive passive galaxies or heavily obscured starbursts. 
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Figure 8. Comparison between typical SEDs for different activity levels. 
All SEDs have been normalized to the K-band flux. Upper panel show the 
GRASIL synthetic SED, while lower panel show the corresponding intrin- 
sic starlight emission. Blue and red lines refer to very active (SSFR > 
10~ 10 yr~ 1 ) and very passive (SSFR < 10 — 14 yr — 1 ) objects in the 
mock catalogue. Cyan, green and magenta lines represents typical SEDs for 
EROs with lO -10 ^" 1 < SSFR < lO -11 ^ -1 , SSFR ~ 10 -11 i/r -1 , 
10 _11 j/r _1 < SSFR < 10 _12 #r _1 respectively. 

To further investigate whether a bimodality is present in the model 
ERO population, we define the color excess E(B — V) as: 

E(B — V) = (B — V) — (B — V)° (1) 

where [B — V)° represents the intrinsic color computed on the 
rest-frame synthetic SEDs before dust extinction, while (B — V) is 
computed after dust extinction. In fig. [5] we show the 2D probabil- 
ity density of control and ERO galaxies in the SSFR - E(B — V) 
color space. Contours are as in figure [4] Also in this space, ac- 
tive and passive EROs do not populate separate regions, but they 
define a unimodal distribution, very similar to that of the overall 
galaxy population. Passive galaxies are characterized by low atten- 
uations as expected, and color excesses range up to a maximum 
of ~ 0.4 for active galaxies. EROs do not populate the region of 
strongly star-forming galaxies with E(B — V) < 0.2. We find 
a clear bimodality when we consider the distribution of the frac- 
tion of galaxies with very red colors, which shows two broad peaks 
for passive galaxies with E(B — V) < 0.2 and very active ones, 
SSFR > 10 _10 2/r _1 and E(B - V) > 0.2. This result does 
not give suppo rt to the assumption of high attenuations for active 
EROs made bv lNagamine et all d2005l) . In fig.[6]we show the loca- 
tion of control and ERO galaxies in the space defined by the mass- 
weighted ages of the stellar populations, normalized to the age of 
the Universe at the observation redshift, and E(B — V). EROs 
tend to avoid the region at low age and attenuation, and again their 
fraction show a bimodal trend in E(B — V), but no bimodality is 
visible in the overall distribution. Finally, we show in figure [7] the 
location of EROs in the stellar mass-metallicity relation. EROs tend 
to selectively populate the high metallicity region. This is the result 
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of the combination of two effects, as can be seen from the smaller 
panels. Passive galaxies populate only the high metallicity region, 
as expected from evolved stellar populations. The control popula- 
tion of active galaxies shows a larger spread, but only sources with 
high metallicity reach the high levels of obscuration required to be 
selected as an ERO. 

The distinction between active and passive EROs is an ob- 
servational crucial point. From figure [4] we see that SSFR values 
for EROs have a broad and unimodal distribution, so that moving 
the cut around the assumed value oflCP 11 yr _1 will not change 
much the relative fractions of the two classes. To better illustrate 
the continuity of EROs along the SSFR parameter, in fig. [8] we 
show five typical SEDs (in the observed frame) for EROs in our 
mock catalogue with various SSFR values. SEDs are normalized 
at the A'-band wavelength, in order to highlight the differences in 
SED shapes. We choose two extreme examples of active and pas- 
sive galaxies (SSFR > 10" 10 j/r _1 and SSFR < 10~ 14 2/r _1 ) 
and three objects with intermediate SSFR; in particular the green 
line represents an object with SSFR close to the activity thresh- 
old, while the cyan and magenta lines refer to objects with activity 
respectively higher and lower than the threshold by less than a fac- 
tor of 10. A and J central wavelengths are shown by dotted and 
dashed lines, while at central A'-band wavelength we report the up- 
per limit for ERO selection. This comparison highlights the contin- 
ually varying properties of SEDs along the SSFR sequence, with no 
clear difference in the predicted spectral properties in the Near-IR. 
Larger differences are present at longer wavelengths (A > 10/im), 
so observations in the Mid-IR and, especially, in the Far-IR can 
provide a clear distinction between active and passive EROs. 

Such observations can in principle be carried out by the Her- 
schel Far-IR space telescope. We show in fig.|9]the predicted rela- 
tion between the SSFR of EROs and the ratio between A'-band and 
monochromatic fluxes in four bandpasses covered by the instru- 
ments on board Herschel, in our smaller mock catalogue where IR 
fluxes are computed using full GRASIL calculations. All four colors 
show a clear correlation with SSFR; we find much weaker trends 
for fluxes at shorter wavelengths, as those currently probed by the 
Spitzer space telescope. These correlations are stronger at high lev- 
els of specific star formation, while passive galaxies split into two 
distinct sequences. This splitting relates to the modeling of star 
formation in bulges or discs: the upper sequence is dominated by 
galaxie s with bulge-to-total mass ratios > 0.6 (as in lFontanot et al.l 
l2009al) . denoted by red crosses while cyan squares correspond to 
disc-dominated galaxies. We checked that this splitting is mainly 
due to the presence of a population of passive and high-metallicity 
discs, that give rise to the lower branch. Getting back to active 
galaxies, it is evident that a Far-IR follow-up of the already-known 
EROs would be of great importance in order to have a better under- 
standing of the properties of these objects, and in particular on their 
star formation rate. We present in fig. [10] the predicted integrated 
number counts for passive and active EROs (red dashed and blue 
dot-dashed lines respectively). In each panel, we also mark the sur- 
face densities corresponding to one detection in the Herschel Key 
Programme^ (GOODS-Herschel, HERMES, ATLAS, see caption 
for symbol coding). The result is not exciting: we predict that it 
is unlikely that these surveys will detect a large number of EROs. 
Incidentally, in the same figure we report our prediction for the in- 
tegrated number counts of the whole galaxy population (thick solid 
line). Our differential number counts are in good agreement (dif- 
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Figure 9. Predicted relation between the SSFR and the the ratio between 
the K-band and the monochromatic fluxes in four Herschel passbands, nor- 
malized to the A'-band flux, for the EROs sample. Red crosses and cyan 
squares correspond to bulge- and disc- dominated model galaxies (splitted 
on the basis of their bulge-to-total ratio, with a threshold value of 0.6). Dot- 
dashed line represents the assumed threshold between active and passive 
EROs. 



ference below 0.2 dex) with the results of lLacev et alj d2009l) for 
infrared fluxes S > 1 mJy; at fainter fluxes MORGANA predicts a 
higher abundance of sources. It will be soon possible to compare 
this prediction with the forthcoming Herschel observations. 

3.2 Evolutionary sequences in morgana. 

In this section we focus on the relation between EROs (defined as 
objects with A — A' > 5.0 and A < 22), sub-mm galaxies (defined 
as galaxies with a predicted 850 ^im flux /sso^m > 0.5 mJy) and 
local massive galaxies (M» > 1O 11 M at z = 0). 

In order to use IR templates in place of the more expensive 
GRASIL, we define sub-mm bright galaxies as those objects with 
/850/j™. > 0-5 m Jy> a uux where number counts based on the two 
methods are very similar (sec. I2.lt and only slightly higher that 
the limit reached in deep fields with currently available sub-mm 
telescopes. 

For these three categories we follow their evolution in the red- 
shift range 0.5 < z < 4.5, with a time sampling of 100 Myr; for 
simplicity we restrict to the main progenitor branch of the galaxy 
merger tree. We check to what extent and in what time order a 
galaxy belonging to one category is observed as a member of an- 
other category. 

As a first figure, the fraction of EROs that are seen at the same 
time as sub-mm sources is low, 0. 1 per cent, while the fraction of 
sub-mm galaxies that are seen as (active) EROs is 5 per cent; then, 
the overlap between the two categories is rather small. In fig.ll llwe 
show some diagnostics of the evolutionary path of galaxies in some 
specific category as a function of their galaxy stellar mass. The 
vertical arrows mark the mass limits within which the diagnostic 
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Figure 10. Predicted number counts in four Herschel passbands. The thick 
solid, red dashed and blue dot-dashed lines refer to the whole popula- 
tion, the passive and active EROs respectively. The symbols mark the sur- 
face density corresponding to the detection of one object in the GOODS- 
Herschel (triangle), PEP GOODS-S (square), HERMES Level-1 (star), Her- 
mes Level-2 (asterisk) and ATLAS (circle) surveys 

has sufficient statistics (only mass bins with at least 10 objects are 
considered). Panel (a) of Fig. II II shows the fraction of EROs that 
are classified as sub-mm galaxy at an earlier epoch as a function 
of their stellar mass at the ERO epoch; the three histograms give 
this quantity for active, passive and all EROs (as specified in the 
legend). Almost all EROs with stellar masses > 10 11 M have 
had a sub-mm ancestor as a main progenitor, irrespective of their 
level of activity. The connection is soon lost at lower masses. In 
panel (b) of fig. [TT] we select z > 1 sub-mm sources and check 
whether they are predicted to be the main progenitors of a (passive 
or active or either) ERO. Here the stellar mass refers to the sub- 
mm selected galaxies. The resulting fraction is high and does not 
depend on stellar mass forM* > lO 9 75 M .This means that most 
sub-mm galaxies have an ERO descendant, which is in most cases 
an active ERO, while only roughly half of them have a passive ERO 
descendant. 

We then consider the expected relation of the two high-redshift 
categories with local massive galaxies. Panel (c) of fi sure [TTI shows 
the fraction of EROs and sub-mm galaxies that are main progeni- 
tors of a z = massive galaxy, as a function of their stellar mass at 
the observation redshift. More massive EROs and sub-mm galaxies 
are more likely to be main progenitors, and the relation is tighter for 
EROs, while sub-mm galaxies have a relatively modest, ~ 30 per 
cent, probability. A different result is obtained when the main pro- 
genitors of local massive galaxies are considered, as in panel (d) of 
figureQj] here we report the fraction of local massive galaxies that 
have an ERO or sub-mm main progenitor, as a function of z — 
stellar mass (obviously larger than 10 1 M©). The main progenitor 
of a massive galaxy used to be an ERO more than 80 per cent of 
times and a sub-mm galaxy ~ 70 per cent of times; as a matter of 
fact, each object may be seen as an ERO or sub-mm galaxy at many 



redshifts. So, the low probabilities shown in panel (c) imply that 
many massive EROs or sub-mm galaxies end in massive galaxies 
not as main progenitors. In other words, consistently with the hier- 
archical clustering scenario, many massive galaxies will have not 
one but several ERO and/or sub-mm progenitors. However, panel 
(d) shows also that only ~ 20 per cent of the main progenitors of 
local massive galaxies were both sub-mm galaxies and EROs (in 
this order) in their past. 

The results in fig.QT]clearly depend on the adopted flux limit 
of sub-mm bright galaxies. Lowering (or raising) this flux has the 
effect of decreasing (increasing) the mass cutoff in panel (a), thus 
increasing (decreasing) the fraction of all EROs with sub-mm pro- 
genitors at intermediate masses. In panel (d), the fractions of sub- 
mm and ERO+sub-mm main progenitor increase (decrease) ac- 
cordingly. Instead, the fractions shown in panels (b) and (c) are 
hardly affected as long as the flux limit is below 2 mJy. 

In fig. [12] we consider the evolution of 5 representative model 
objects, chosen between the ~ 25 per cent explicitly following the 
sub-mm - EROS - massive galaxy evolutionary sequence. In the 
four panels, from top to bottom, we show the total stellar mass 
(symbols mark the 2 = mass), the specific star formation rate, 
the R — K color and the 850/im flux. The diversity of evolutionary 
tracks is evident from the figure, as well as the complexity of star 
formation histories. In most cases the two phases are well separated 
along the redshift evolution of the object. Some model galaxy un- 
dergoes repeated independent ERO (and sub-mm) episodes along 
their evolution; some of them can be seen both as an active and a 
passive ERO at different redshifts. Even if in most cases the ERO 
stage follows the sub-mm bright period, it is also possible that, as 
a result of late star formation activity, a model galaxy may be rec- 
ognize as a sub-mm bright galaxy, even if it was a passive ERO 
at previous cosmic times. Interestingly, active EROs tend to have 
low sub-mm fluxes, consistently with the low fraction of galaxies 
that are at the same time ERO and sub-mm. These results show 
that the galaxy evolution sequences in our galaxy formation model 
based are more com plex compared to the simple one obtained by 
iGranato et all j2004l) . that predict that massive galaxies go through 
sub-mm, quasar, ERO and passive phases. 



4 CONCLUSIONS 

In this paper we have discussed the statistical properties of the ERO 
population as predicted by the MORGANA galaxy formation model, 
in terms of number counts, redshift distributions and active frac- 
tions of R — K > 5 sources. We generated mock photometric cat- 
alogues by interfacing MORGANA with the spectro-photometric + 
radiative transfer code GRASIL, and created large samples of model 
galaxies (so as to have good statistics on the rare EROs). We used 
both the Padova and the M05 SSP libraries; the latters include the 
contribution of TP-AGB stars. We combined the resulting UV-to- 
Near-IR synthetic SEDs, computed w ithout the expensive Far-IR 
emission, with the iRieke et all d2009t) IR template library to esti- 
mate IR fluxes longward 3fim. We then selected samples of EROs 
(and sub-mm galaxies) by applying the corresponding color and 
flux cuts. 

We found that our standard model with Padova SSP library is 
able to reproduce the overall number counts and redshift distribu- 
tion of the R — K > 5 sources; the agreement worsens when red- 
der cuts are considered. Other discrepancies, like the excess of faint 
EROs at deep K fluxes and the underestimate of the number of dis- 
tant passive EROs, are consequences of well-known points of ten- 
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Figure 11. Panel (a): Fraction of EROs with a sub-mm main progenitor with fssoum > 0-5 m Jy m MORGANA. Thick solid, red dashed (diagonal texture) and 
blue dot-dashed (cyan shading) histograms refer to the total and passive and active EROs respectively. Panel (b): Fraction of sub-mm galaxies that are the main 
progenitor of an ERO in MORGANA. Thick solid histogram refers to the total number of sub-mm galaxies, while the red dashed (diagonal texture) and blue 
dot-dashed (cyan shading) histograms indicate the properties of the descendant as indicated in the caption. Panel (c): Fraction of EROs (solid histogram and 
yellow shading) and sub-mm galaxies (green dot-dashed histogram and diagonal texture) that are main progenitors of a local massive galaxy (M, > 10 11 Mq) 
as a function of their stellar mass. Panel (d): Fraction of massive (M* > 10 iW©) galaxies at z < 0.2 with a EROs (solid histogram yellow shaded area), 
a sub-mm progenitor (/gso^m > 0.5 mJy green dot-dashed histogram and diagonal texture) or both (thick solid histogram). In all panels the vertical arrows 
indicate the mass limits for the reference sample of EROs, sub-mm and massive galaxies. 



sion of the model with data, like the dearth of ma s sive z > 2 galax- 
ies jFontanot et all |2007| ; Icirasuolo et alj|2008l : iMarchesini et alJ 
200^7 or the excess of z ~ 1 small galaxies l IFontanot et al.l 
2009bh : this analysis does not reveal new discrepancies. Interest- 



ingly, MORGANA is able to roughly reproduc e the substantial con- 
tribution (of the ord er of ~ 50%, see e.g. ICimatti et alj|2002al ; 
iMivazaki etai]l2003l) of active galaxies to the global ERO popula- 
tion; here we define active galaxies as those with SSFR > 1CP 11 
Gyr. We ascribe this suc cess to the higher SSFR predicted by th e 
model at higher redshifts dFontanot et alj200 7: Santini et alj2009h . 

We showed that the ERO number counts increase consider- 
ably when using the M05 SSP library. The addition of TB-AGBs 
gives a strong boost to the restframe SEDs of active galaxies long- 



ward 6000 A, reddening them and thus increasing the number of 
active EROs. This leads to a good level of agreement even at very 
red cuts (R — K > 6), though it does not completely solve the 
lack of very red passive galaxies at higher redshift. This highlights 
the importance of an accurate modeling of galactic synthetic SEDs 
when comparing the prediction of theoretical models of galaxy 
formation to observations, and confirms that TP-AGBs are a rel- 
evant ingr edient for the correct reproduction of the high-redshift 
Universe dMaraston et~ai1l2006|) . The details of the modeling of 
dust attenuation also play a role: we checked that both increas- 
ing the fr action of cold gas l ocked into molecular clouds (as sug- 
gested bv lLo Faro et al.l2009l to reproduce the properties of Lyman 
break Galaxies at z > 4) or using the analytical approach pro- 
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Redshift 

Figure 12. The redshift evolution of the total mass (panel a), SSFR (panel 
b), observer frame R-K color (Panel c), and observer frame 850/^m flux 
(Panel d) for 5 representative galaxies. In each panel the same greyscale 
corresponds to the same object. Empty symbols correspond to the 2 = 
mass of the 5 galaxies. 



posed bv lDe Lucia & Blaizoj d2007h based on an universal attenu- 
ation law slightly reduces the number of active EROs ( £ 0.3 dex 
at K ~ 20). 

We then considered the physical properties of the predicted 
ERO population with respect to a control sample of non color- 
selected galaxies with A' < 22 and 1 < z < 2. Overall the physical 
properties of EROs do not differ significantly from those of a typi- 
cal galaxy in the same redshift range, apart that active galaxies with 
low dust attenuation and low metallicities are excluded. We do not 
see any signature for a significant bimodality in the ERO popula- 
tion: the probability of a galaxy to be selected peaks in two regions 
of the SSFR-E(B — V) space, corresponding to moderately atten- 
uated passive galaxies and strongly attenuated (E(B — V) > 0.2) 
active galaxies, but this trend is not visibl e in the overall ERO 
distrib ution. This result is in agreement with lGonzalez-Perez et al.l 
d2009h . Deep fields in the Far-IR, from 100 to 500 fim, are the best 
way to constrain the SSFR of these objects, but presently planned 
surveys with Herschel will detect a small number of EROs at best. 

We also tested the possible relation between EROs and both 
the sub-mm (fssopm > 0.5 mJy) galaxies and the local massive 
(M, > 10 11 Mq) galaxies. This was done to test the suggestion 
that these different populations reflect well-defined stages of a typ- 
ical evolutionary sequence of massive galaxies, like sub-mm galax- 
ies evolving into EROs and later into massive passive galaxies. We 
found that most massive EROs have a sub-mm bright main progeni- 
tor and evolve into a massive galaxy at z = 0, and that many z — 
massive galaxies have an ERO or a sub-mm main progenitor, but 
only a minor fraction, ~ 25 per cent, follow the sub-mm - ERO - 
massive galaxy sequence. Moreover, massive galaxies may easily 
have more than one sub-mm and/or ERO progenitor. In fact, star 
formation histories predicted by MORGANA are far more complex 



and exhibit a great degree of diversity, with only a small sample 
following the well defined sequence described above. 

Our work shows that EROs can be used as a powerful con- 
straint for theories of galaxy formation and evolution, but a better 
and deeper understanding of the distribution of EROs in a configu- 
ration space defined by the physical quantities, like star formation 
activity, stellar mass, metallicity, stellar population ages and dust 
obscuration, is needed. Unfortunately, planned surveys with Her- 
schel may not reach the combination of depth and area required for 
detecting and analysing statistically significant samples of EROs. 
On the other hand, the same instruments may provide the required 
spectroscopic and photometric follow-up in order to better under- 
stand the real nature of already detected EROs. 
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